LIPSCHITZ STABILITY OF AN INVERSE BOUNDARY VALUE 
PROBLEM FOR A SCHRODINGER TYPE EQUATION 

ELENA BERETTA*, MAARTEN V. DE HOOPt, AND LINGYUN QIU* 

Abstract. In this paper we study the inverse boundary value problem of determining the 
. . . potential in the Schrodinger equation from the knowledge of the Dirichlet-to-Neumann map, which 

^^ is commonly accepted as an ill-posed problem in the sense that, under general settings, the optimal 

stability estimate is of logarithmic type. In this work, a Lipschitz type stability is established 
f^^ 1 assuming a priori that the potential is piecewise constant with a bounded known number of unknown 

^vj ' values. 

> 

Q ■ 1. Introduction. In this paper, we investigate the stabihty for the inverse 

boundary value problem of a Schrodinger equation with complex potential, q{x) say. 
This encompasses the Helmholtz equation with attenuation, when q{x) — cj^c~^(x), 
where c denotes the speed of propagation and w is the frequency, which can be com- 
plex. In fact, the imaginary part of U!c''^{x) characterizes the attenuation of waves in 
the medium. 

Ph ■ We begin with formulating the direct problem. Let u e H^{n) be the weak 

■^^ I solution to the boundary value problem, 

'5 ■ n 1^ r i-A + q{x))u^ 0, xen, 

c^ ; ^'-^> \ u^ g, xedQ, 

where fl C M", rt > 2 is a bounded connected domain, q S L°°(ri) is a complex- valued 
function and g is prescribed in the trace space H^^^{fl). The Dirichlet-to-Neumann 
^l^J ■ map is the operator A, : iJ^/^(Sl) — > H~^/'^{Q) given by 
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where i^ is the exterior unit normal vector to dfl. 

The inverse problem that we consider, consists in determining q when Ag is known. 
This problem arises in geophysics, for example, in reflection seismology assuming a 
description in terms of time-harmonic scalar waves. The topic of this paper is the 
issue of continuous dependence of q from the Dirichlet-to-Neumann map A^. The 
continuous dependence is of fundamental importance for the robustness of any re- 
construction, as well as for the development of convergent iterative reconstruction 
procedures starting not too far from the solution (cf. [5]). More precisely, it has been 
proved that Landweber iteration reconstruction methods converge if the continuous 
dependence for the inverse problem is of Holder or Lipschitz type. 

From the work of |10| . it is evident that for arbitrary potentials q, Lipschitz stabil- 
ity cannot hold. Motivated by, and following analogous results in electrical impedance 
tomography (KIT, cf. [3l|4]), here we study conditional stability when a-priori infor- 
mation on q is assumed. We consider models with discontinuous potentials to ac- 
commodate realistic reflectors. Specifically, we consider the space spanned by linear 
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combinations of N characteristic functions. More precisely we consider potentials of 
the form 

N 

J=l 

where qj, j = 1, . . . N are miknown complex numbers and Dj are known open Lipschitz 
sets in R". Moreover, we consider the case of partial boundary data, that is, we can 
restrict the collection of measurements to only a part of the boundary. We refer 
to ^13j for a review of recent uniqueness results. Here, we prove Lipschitz stability 
with a uniform constant, which depends on N and on the other a-priori parameters 
of the problem. We will show that the Lipschitz constant grows exponentially with 
the dimension, iV, of the space of potentials. The method of proof follows the ideas 
introduced in Alessandrini and Vessella and relies on quantitative estimates of unique 
continuation of solutions to elliptic systems and on the use of singular solutions and 
of their asymptotic behaviour near the discontinuity interfaces. Compared to the case 
of the real or complex conductivity equation in the case of the Schrodinger equation 
we are able to derive our result relaxing the assumptions of regularity on dDj that are 
assumed to be Lipschitz. Furthermore, taking advantage of the regularity of solutions 
and of its gradient inside the domain ft we find a better dependence of the stability 
constant on N. 

The outline of the paper is as follows. In the next section we state all the as- 
sumptions and the main result. In Section [3l we give a summary of known regularity 
results connected to Schrodinger equation with complex potential, and some prepara- 
tory lemmas concerning the existence and asymptotic behaviour of singular solutions. 
Section |4] contains the proof of our main theorem. We first show the proof for n = 3 
and then modify it to the other cases. For the structure of the main proof we charac- 
terize the rate of blow-up of the singular functions finding lower and upper bounds in 
terms of the distance of the singularity from the interface of the subdomains. More 
precisely, to derive our main result we first establish that the singular function sat- 
isfies a lower bound in terms of the distance of the singularity from the interface. 
Secondly, by using quantitative estimates of propagation of smallness we derive also 
an upper bound for the singular function. Last but not least, we make use the value of 
a bounded non-decreasing function at some particular point to prove that either the 
result of the main theorem can be deduced directly or a recursive inequality (|4.22p 
must hold true. The recursive inequality also leads to the desired result. In Sec- 
tion [5] we demonstrate by an example that the Lipschitz constant grows exponentially 
with the dimension of the space of potentials. This example is constructed from its 
analogue in electrical impedance tomography [IV . 

2. Main result. 

2.1. Notation and definitions. We denote by n the space dimension. For 

every x G M", we set x = (a;',a;„) where x' € M"-^ for n > 2. With Br{x), B'j^{x') 
and Qr{x) we denote the open ball in M" centered at x of radius R, the ball in R"~^ 
centered at x' of radius i?, and the cylinder B'j^{x') x (a;„ — R,Xn + R), respectively. 
For simplicity of notation, Bfj{0), -BJ^(O) and Qb.{0) are denoted by Bfj, B'j^ and Qr. 

Definition 2.1. Let fl be a bounded domain in R". We say that a portion E 
of dfl is of Lipschitz class with constants ro,L > if, for any P G E, there exists a 
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rigid transformation of coordinates such that P ~ and 

nnQra ^ {(a;', Xn) G Qrg I Xn > H^')} 

where (j) is a Lipschitz continuous function on B'^ with (/)(0) — and 

We shall say that fl is of Lipschitz class with constants rg and L, if dfl is of Lipschitz 
class with the same constants. 

Definition 2.2. Let fl be a bounded open subset o/ M" and of Lipschitz class 
and Yi be a open portion of dO.. We define Hco (S) as 

Hli^T.) = {g e H^'^{d^) I SUPP5 C S} 

and Hco (S) as the topological dual of Hd (S); we denote by (•,•) the dual pairing 
between Hco (S) and Hco (S). 

Definition 2.3. Let fl be a bounded open subset ofW^ and of Lipschitz class, 
T, be a open portion of dfl and q € L°°(il). Assume that is not an eigenvalue of 
(—A + q) with Dirichlet boundary conditions in fi, i.e., 

{u G Hl{^) I (-A + q)u = 0} =: {0}. 
For any g G Hco (S), let u G H^{Vl) be the weak solution to the Dirichlet problem 

r {-/^ + q{x))u^ 0, xGfi, 
^ > \ u^ g, xedfl. 

We define the local Dirichlet-to- Neumann map Aq as 

du 
where v is the exterior unit normal vector to dQ. 



With il being a bounded open set, with C"'^ boundary, the set of the eigenvalues 
of (—A + q) with Dirichlet boundary conditions is a discrete subset of C, and hence 
can be avoided. 

We observe that A, can be identified with the sesquilinear form on Hco (S) x 
Hli^{Y), defined by 

(Af)5,/)= f {^u-Vv + quv)dx, yf,geHU\j:), 



where u is the solution to (j2.ip and v is any function in H^{Q.) such that i; \dn= f ■ 
This definition is independent of the choice of v. Let v\, V2 be two different functions 
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in H^{Q,) such that vi \gn— ^2 \dn— f ■ Then, since w — Vi — V2 ^ H}^{Q?), and m is a 
solution, we have 

/ (Vu ■ Vwi + 9UU1) dx — / {V u ■ 'SI V2 + quv2) dx 
Jn Jn 

~ / (VuVw + quHj) dx = 0, 
Jn 

using integration by parts. We denote by || • |L,rri/2,^N f^-l/2,^^^ the norm defined as 

2.2. Main assumptions. Our assumptions on il and q{x) are 
Assumption 2.4. fi C K" is a bounded domain satisfying 

\n\ <A 

Here and in the sequel \fl\ denotes the Lehesgue measure offl. We assume that dfl is 
of Lipschitz class and we fix an open portion S of dCl which is of Lipschitz class with 
constants tq and L. 

Assumption 2.5. The complex-valued function q{x) satisfies 

\\q\\L^m < B, 
where B is a positive constant, and is of the form 

N 

where qj,j — 1, . . .TV are unknown complex numbers and Dj are known open sets in 
M" which satisfy the following assumption. Moreover, we assume that is not an 
eigenvalue of — (A + g) with Dirichlet boundary conditions in fl. 

Assumption 2.6. The Dj,j ~ 1, . . . , N, are connected and pairwise non- overlapping 
open sets such that UjLiDj ~ fl and dDj are of Lipschitz class. We also assume 
that there exists one set, say Di, such that dDi D dQ contains an open portion Si 
of Lipschitz class with constants ro and L. For every j G {2,...,A^} there exist 
ji, . . . , JM e {1, • • • , A^} such that 

and, for every k ~ 1, . . . , M , 

contains a non-empty open portion S^ of Lipschitz class with constants tq and L such 
that 

SiCS, 
Efe C r^, \Jk^2,...,M. 
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Furthermore, there exists Pk G S^,, at which -D^-i satisfies the interior ball con- 
dition with radius ^, and a rigid transformation of coordinates such that P^ = 
and 

Sfe n Qro/3 = [x € Qra/3 I Xn = <Pk{x')}, 

Djk n Qro/3 = {a; e Qra/3 I Xn > (f>k{x')}, 

Djk-l n Qrn/3 = {a; e Qro/3 | Xn < (l>k{x')}, 

where (j)k is a C"'^ function on B'^ ,3 satisfying 



0fc(O)=O 



and 



ll<^fe||c».i(B;^/3) <i- 

for simplicity, we call Dj-^ , ■ • • , ^jAf '^ chain of domains connecting Di to Dj. 















Dk 








Pk 




I.. 




Pk-i 





In the further analysis, for simphcity of notation, we also use the constant ri — j^. 

2.3. Statement of the main result. The main result of this paper is stated 
as follows. 



Theorem 2.7. Let f2 satisfy Assumption \2.4\ and q^''\k — 1,2 he two complex 
piecewise constant functions of the form 



N 



,(fe) 



(-) = E 



jfc). 






which satisfy Assumption \2.5\ and Dj,j — 1,...,N satisfy Assumption \2.6i Then, 
there exists a constant C — C{n, tq, L, A, B, N), such that 



(2.2) 



lk« -g(^'llL~(o) < GIIAf' - Af)||^(^v.(^^^^-v.(^)), 



where kf^ = A^f,^, fork = 1,2. 



3. Preliminary results. In this section, we state some results which will be 
used in the proof of our main stability result. 

Proposition 3.1. Let il be a bounded Lipschitz domain in W\ q e L^{fl) 
complex valued potential , f £ LP{ri) and g £ W ^p'^{dQ) with 1 < p < oo. Assume 
that is not a Dirichlet eigenvalue for the operator —A + q in fl. Then there exists 
a unique solution u G W'^''P{Vt) to the problem 

.g^N f {-A + q{x))u^ /, xef], 

^ ' ' \ u — g, X e 951, 

Moreover, 

(3.2) \\u\\w^..(^n) < C (ll3ll^.-i,.(^^) + Wfhnn)^ 
where C depends on n,Q and \\q\\L^(n)- 

The proof is a consequence of the of existence of a VF^^^(ri) function w such that 
w = g on dn and such that ||w||i^2,p(q) < C\\g\\ 2-i,p and of the Fredholm 

alternative; see for example Theorem 3.5.8 in Feldman and Uhlmann's notes [7]). For 
reader's convenience, we also note the following Proposition 13 . 21 without proof, which 
we use for the low dimension cases. 

Proposition 3.2. Let Q, be a bounded Lipschitz domain in R", q e L°°{fl) 
complex valued potential , f G H~^{^) and g £ H^''^{dil). Assume that is not a 
Dirichlet eigenvalue for the operator — A + g in fl. Then there exists a unique solution 
u £ H^{fl) to the equation iS.l]) . Moreover, 

(3.3) \\u\\HHn) < C (MHU^on) + 1I/1Ih-(o)) 
where C depends on n,n and \\q\\L°°in)- 

Our approach follows the one of Bcrctta and Francini[4 , which is for the FIT 
problem with complex conductivity, of constructing singular solutions and of study- 
ing their asymptotic behavior when the singularity approaches the interfaces Sfc. This 
method was originally introduced by Alessandrini and Vessella in the real- valued con- 
ductivity case [3]. To construct singular solutions for the FIT problems, the Green's 
function plays a crucial role. In our case, we also use the Green's function to treat 
the case of high dimension {n > 4) and a first order derivative of Green's function 
needs to be used for lower dimension (n — 2,3). In the following propositions, we 
discuss the existence and behavior of the Green's functions {n > 4) and a first order 
derivative of the Green's function (n = 2, 3) when q satisfies Assumption 12.51 We are 
especially interested in their asymptotic behavior near the C°'^ interface Efc. 

Before doing this, we need to extend our original domain. We consider Si and 
recall that up to a rigid transformation of coordinates we can assume that Pi = and 

{W"\n) n Bro = {{x', X„) e Br„ I Xr, < (t>{x')} 

where (/> is a Lipschitz function such that 0(0) — and ||0||c".i(S' ) ^ L. Then we 
extend 51 to J7o = 51 U D^ by adding an open set Dq defined as 



Do^ixc (M"\51)nBr 



12. 
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< -^"Tq, \xi\ < -ro, « = l,...,n- 1 
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It turns out that Qq is of Lipschitz class with constants ^ and Li, where Li depends 
on L only. We define 



Ko^{xeDo |dist(x,I]i)> y} 



with dist(iiro, dil) > ^. We extend q{x) defined on il by setting it equal to 1 in Dq. 
For simplicity of notation we still denote this extension by q{x). 

We consider any subdomain in 51 and the chain of domains connecting it to Di. 
For simplicity let us rearrange the indices of subdomains so that this chain corresponds 
to Dq, Di, . . . , Dm, M < N . Let S — U^LqDj and K he & connected subset of S with 
Lipschitz boundary such that K D dDj = Ej U Sj+i for j = 1, 2, . . . , M, Ko C K and 
dist(/^,a5\{SA,+iUEi})>f|. 
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In the following, we shall use C to denote positive constants. The value of the con- 
stants may change from line to line, but we shall specify their dependence everywhere 
where they appear. For n > 4, let F denote the fundamental solution associated with 
the Laplace operator. In the proof of Theorem 12. 7[ we will need to estimate G — F 
from above in terms of variable-interface distance r to a power, which is smaller than 
the order of the singularity of F. Since, for high dimension cases(n > 6), F(-, y) does 
not belong to 7J~^(ri), we need to employ L^ estimate of the solutions here. Note 
that F(-,y) belongs to LP(fi) for any 1 < p < ^. 



Proposition 3.3. Let the complex-valued function q G L°°(J7o) satisfy Assump- 
tion \2.5\ and n > A. For y G flo, there exists a unique function G{-,y) continuous in 
flo\{y} such that 



(3.4) 



VG(-, 2/)V0 + qG{; y)4> = 0(y), V0 e C^{n). 



Furthermore, we have that G{x, y) is symmetric, that is, 
(3.5) G{x, y) ^ G{y, x), x,y e ilo, 

and the following estimates 



(3.6) 



\\Gi-,y)\\L^no\BAv)) ^ C\lnr\ 
\\G{;y)\\LHno\BAy)) < <^^'"* 



r < -dist(y,5rJo), n = 4 



r < -dist{y,dn„), n>5 



and 



C , 4 < n < 7, 



n = 8, 



(3.7) l|G(-,2/)-r(.,y)|U.(f,„)<<( |ln(dist(2;,U,ti9i?,))| , 

dist{y,uf^^dDjf'^ , n>9, 

for dist(?;, dflg) > ^ , hold true, where the constant C depends on the constant in 
Proposition \3.1\ 

Proof. Assume that y belongs to some sub-domain Z?,„ which q equals to a 
complex constant qm inside. Let H{x,y) denote the outgoing fundamental solution 
of Helmholtz equation 

{^A + q,n)H{x,y)=S{x,y), a; £ R", 

i.e., 



(3.8) Hix,y) 



qh 



(n-2)/4„(l) (a^''^\x-y\\ 



(n-2)/2 



4i(27r)("-2)/2|2._y|(n-2)/2 ' 



where iJ„ denotes Hankel function of the first kind. We consider G(x, y) = H{x, y) + 
u](x,y), where to solves 

/g gN f (-A + q)uj = {q„i - q)H, in f^o, 

^ ' ' [^ cj = —H, on OHq. 

Note that qm~q vanishes in Dm- Hence {qm — q)H belongs to L°°{^q). By using the 
asymptotic behavior of the Hankel function near the origin [T2] , we obtain that 



|(g„ -q)H{x,y)\ < ^ ^^ 



, \x-y\< dist{y,dDm), 

Clx-yl^-" , \x^y\>distiy,dDra), 



for some positive constant C. We observe that the order of the singularity of a;(a:, y) is 
always lower then the fundamental solution H{x, y). To be more precise, by applying 
Proposition 13.11 with p = ^^ and Sobolev embedding theorem, we conclude that 

(3.10) M;y)\\mn,) <CM;y)\\ ^^^ 

r C , 4 < n < 7, 

<C||(g„-g)i/(-,2/)|L ^^^^ < J |ln(dist(y,9An))| , n = 8, 

ydist{y,dDm) ^ , n>9. 

Then, using the asymptotic behavior of the Hankel function again and the inequality 

\\G\\L^(n„\B^(y)) < \\^\\L^no\Br(y)) + \\H\\mno\BAy)) 



we immediately get p.6p . 

Let r(-) stand for the Gamma function. Noting that 

Hi;y),ri;y)eC°^{no~{y}) 



and 

H{x,y) -T{x,y) 



n \ 2 j 4i7r("-2)/2l'^ 2/1 n(2 - n)7r«/2 1^ ^1 

= 0, 

as |a; — j/| goes to 0, we conclude that |r(-,j/) — 7J(-,y)|is uniformly bounded for all 
y such that dist(y,ar2o) > %■ Then ^^ follows. D 

In both Beretta & Francini's proof [4] and Alessandrini & Vessella's proof [3], the 
blow-up property of a singular function, 

VGi(2/,a;)VG2(a;,y)da;, 

where Uk = ^\ Uj=i Dj and Gi, G2 are functions defined by p.4p for potentials g(^\ 
q^^', respectively, when y approaches the interfaces, is essential. However, in the case 
of the Schrodinger equation, this does not happen if n = 2, 3. Therefore, for n = 2, 3, 
we will introduce a derivative in the point source. For n = 3, let 

r(x,,)- ^^-^^ 



A'k\x — j/P ' 
which is the solution to the equation 

(3.11) ^ ^T{x,y) = ^5y{x) 



dx. 
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Proposition 3.4. Let n = 3 and q G L°°(fio)- For y 6 f^o, there exists a unique 
function G{-,y) continuous in fioMj/} such that 

(3.12) / VG(-,2/)-V0 + (zG(-,y)0=^— 0(y), V(/. S Go°°(f]). 

Furthermore, we have that G{x,y) is symmetric, i.e., 

(3.13) Gix,y) = Giy,x), a;,yGf^o, 
and the following estimates 

(3.14) \\G{;y)\\mno\B,.iy)) < Cr-^ r < ^dzst{y,dn) 
and 

(3.15) ||G(.,y) -r(.,y)|U.(^„) < G, dist{y,dno) > ^. 
hold, where the constant C depends on the constant in Proposition 

Proof. Consider G{x, y) = T{x, y) + uj{x, y), where uj solves 



I 



(3.16) <! (-^ + '^)-= '^r i^^o^ 

^ ' ' io = — r, on o\Iq. 



Since T{-,y) G M^5/4,4/3(g^^)^ ^^ ^ ^^/^(f^o) and -r(-,y) e H^^^{dno), by Propo- 
sition [3H1 p.l6p has a unique solution w e iJ^(ilo) and a; = G — F satisfies the 
estimate 

(3.17) ||^(-,2/)||Hi(ao) <C(||r(.,y)||^v.(a^„) + ||<?(-)r(-,2/)||,/-i(ao)) <C, 
when dist(y,9J7o) > fg- Hence, 

(3.18) \\G{-,y) - r(-,y)||i2(n^\B,(j^)) = ||tj(-,y)|U2(Oo\B,(y)) 

< \\^i-,y)\\H^no\B^(y)) < \\^i-,y)\\m{no) < C- 
With the fact that 

(3.19) ||r(-,y)||L=(ao\B.(y))<C'^^^ r < idist(zy,ar!), 
(I3.18P gives the desired estimate 

(3.20) ||G(•,y)|U2(o„\B,(,))<Cr-^ r < idist(y, 917). 
Finally, again by p.l7p we immediately get 

(3.21) ||G(.,2/)-F(.,y)|U.(^„)<C. 
D 

For n = 2, let 

^, ^ 27r(x2 - y2) 

which is the solution to the equation 



(3.22) -^Y(x,y)^^^6y{x) 



d_ 
dx2 



Proposition 3.5. Let n — 2 and q e L°°(r2o). For y e r^o, there exists a unique 
function G{-,y) continuous in iloMj/} such that 

(3.23) [ (^s/G{;y)-Vcf> + qG{;y)^=^c^{y), \/cf>eC^{n). 



dXn 



Furthermore, we have that G{x, y) is symmetric, that is, 

(3.24) G{x,y)^G{y,x), x,yeno, 
and the estimates 

(3.25) ||G(.,y)|U.(^^,\s.(^))<qinr|^ r < 1 min ('dist(y,af]o), ^ 

and 

(3.26) ||G(.,y)-F(.,y)|U2(,j„)<C, dist(y, 917o) > ^ 
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hold, where the constant C depends on the constant in Proposition 

We omit the proof here, because it follows from an adaption of the proof of Propo- 
sition 13.41 The symmetry of G follows by standard arguments based on integration 
by parts (see for example [6]). 

In the sequel we will derive estimates of unique continuation in K for solutions 
to our equation. A key ingredient to obtain these estimates is the Three Spheres 
Inequality that we will state below and that was proved by [TJ Theorem 3.1]. The 
next two propositions concern Three Sphere Inequalities for our equation. To prove 
it, one interprets the equation {—A + q)u = for a complex function q{x) as a weakly 
coupled system of equations with Laplacian principal part 

(3.27) -AU + QU = 0, 

where C/ is a vector with components the real and imaginary parts of u, that is, 
y(i) = yiu^ y(2) — g^^ and Q is a two by two tensor with elements the real and 
complex part of the potential q, that is, 9^^^ — Dlq and q^-^^ — 3q. We can also write 
the system in the form 

{ -Au(2)+5(l)^,(2)+ ^(2)^(1) = 0. 

In [H Theorem 3.1] the authors prove the validity of the Three spheres inequality for 
elliptic systems with Laplacian principal part. In particular it applies to solutions U 
of (|3.27p and hence also to solutions of {—A + q)u = 0. 

Proposition 3.6. Let u be a solution to the equation 

{-A + q)u^O in Br. 

Then, for every pi,p2,P3, with < pi < P2 < Ps < R, 

(3.28) Ml^b,.^) < Q2||«||2^(B,jlklli^(\3), 

where a = . ^ 6 (0, 1) and Q2 > 1 depends on \\q\\L'=^{BR) , — cind —. 

PI P F 

Remark 3.7. In 11^ Theorem 3.1] the authors prove the validity of the three- 
spheres inequality for elliptic systems with some limitations on the radii. The deriva- 
tion of the inequality for arbitrary radii follows by applying the argument of the proof 
of in Theorem 5.1] choosing Bro{xo) = B^ , G = -8^2 and 17 = B^^ ■ 

Also, we have 

Corollary 3.8. Let u be a solution to the equation 

{-A + q)u = {) in Br. 

Then, for every pi, P2, Ps, with < pi < P2 < P3 < R, 

(3.29) \\u\\l^(b,,) < Q^\M1<^(b,,)\M]7Ab,,)' 
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where (3 = ^'^^tl'' G (0, 1) and Qoo > 1 depends on \\q\\L^(^BR), j[ and ^. 

Proof. We use the local boundedness estimate for u'^^ and u^^\ weak solutfons 
of elliptic equations (see for instance Theorem 8.17]), to obtain that there exists a 
constant C, which only depends on n and HfzllLoorBH): such that 

(3-30) ll"ll^-(B..)< (^3„^^)„/2 ll"IU-(B.3)- 



Then, by Proposition 13.6 



C 



'^^-P2) 



CQ2 ||_^j|Q IL,I|1-Q 

2 P^) 

CQ2 



(3.31) ^ ,p,+p3 „ ,n/2 ll^ll^-(i3p.)ll"llL^(B,3) 

< — ^^ — ^ii?„, r/^iB„j(^'")/^ikii?oo,„ ^lk"l-" 



D 

As a consequence of the Three Spheres Inequality stated in Corollary 13. 8[ we 
derive the following quantitative estimate for unique continuation of solutions to our 
equation. 

Proposition 3.9. Let K and Kq be defined as before, and let v g H^{K) be a 
weak solution to the equation 

{-A + q{x))v = in K. 

Assume that, for given positive numbers Eq, Eq and real number 7, v satisfies 

(3.32) M|l=o(Ko) < £0, 
and 

(3.33) \v{x)\ < (eo + Eq) dist(a;, Y^m+iV , x e K. 
Then the following inequality holds true for every < r < 2ri , 

/ r \ ^-^"' 

(3.34) \v{i)\ < C —^ (eo + E^y^''^^^'' , 

where x = Pm+i ^?'^(-Pa/+i) with v being the exterior unit normal vector to ODm at 

. / 12ri-2r \ 

Pm+1, /? = hi4 ' ''"'' ~ — /S'-^-'-V ^ ('^' ^^ '^"■'^ ^^^ constants Ni and C depend on 
ro, L, A, B and n. 



'"(tft 



Proof. We construct a chain of spheres of radius ri with centers xq^xi, . . . ,Xk 
such that the first is B^^ (a;o) C i?4ri {xq) C Kq, all the spheres are externally tangent, 
and the last one is centered at x^ = Pm+i ~ 3?'ii/(P/vf+i)- We choose this chain so 
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that the spheres of radius 4ri concentric with those of the chain, except the last one, 
are contained in K and have a distance greater than ri away from Y,m+i- Such a 
chain has a finite number of spheres that is smaher than iVi — „ , + 1. 

By CoroUarv 13.81 and (|3.33p . we have 



\\v\\L°-(Br^{xi)) < II«IIl°=(B3,.i(2;o)) 



where C depends on Q^o and ri. By iterated apphcation of Corollary 13.81 to v with 
radii ri, 3ri and 4ri over the chain of spheres, we have, by p.32p . 

where C depends on Qoo and ri. Now, we let x — Pm+i — ri^iPM+i) where r < 2ri. 
Using Corollarv 13.81 again for spheres centered at Xk of radii ri, 3ri — r and 3ri — |, 
we obtain that 

< c 

which completes the proof. D 



{t^T (- + -«)"'-'-. 



Remark 3.10. Let us observe that, in order to apply Proposition \379\ to the 
singular function defined in Section^when n = 2,4, we need to replace the condition 

(EM) by 



(3.35) \vix)\ < (eo + -Bo) |ln(dist(a;,EM+i))|5, x e K. 

By using the same proof technique, we can obtain the same result with ^3.34^ replaced 
by 

(3.36) \y(^i)\ < C ( ^^ ] (eo + £;o)|lnr|- 



4. Proof of the main result. Assume that Dm is the subdomain of the par- 
tition of fl where the maximum of \\q^^' — q^'^'W is realized and let us denote 

(4.1) E = ||gW - <z(^)||l^p,,) = IkW _ q(')U^^ay 

We consider the chain of domains, Dq, Di, . . . , Dm, as before; S, K and Ko are defined 
as in the previous section. We set 



C/o = n,Uk^ n\ U^^^i D,, fc = 1, . . . , M and Wk = ^j=oD 
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Let y (£ K. For dimension n > 4, let Gi{x,y) and G2{x,y) be the Green's function 
related to g'^-* and q'^'^\ respectively, the existence and behavior of which was shown 
in Proposition 13.31 For dimension n = 2,3, let Gi{x,y) and G2{x,y) be a first order 
derivative of the Green's function, the existence and behavior of which was shown in 
Propositions 13.51 and 13.41 respectively. We define 



(4.2) Skiy,z)= {q<^''>-q^'^)ix)Gi{x,y)G2{x,z)dx. 

JUk 

By Proposition 13. 3[ 13.41 and 13. 5[ there exist a constant C such that 
(4.3) 

\Sk{y,z)\<CE\\n{d\sl{y,Uk))Hd\st{z,Uk))\^, y,zeKnWk, n = 2,4; 

\Sk{y,z)\<GE {dist{y,Uk)dist{z,Uk)rK y,zeKnWk, n = 3; 

\Sk{y,z)\<GE {dist{y,Uk) dist{z,Uk))^'^, y,zeKnWk, n>5. 

We focus on n = 3 first; we will discuss the adaptation of the proof for the case 
n = 2, 4 and n > 5 at the end of the proof. 

Lemma 4.1. For every y,z £ K n Wk, we have Sk{-,z), Skiy, •) G H^{K n Wk) 
and 

(4.4) (-A + g«)5fc(.,z) = 0, {-A + q^^'>)Sk{y,-)^0 mKDWk. 



The proof of this Lemma follows from the symmetry of Gi (i = 1, 2) and changing 
the order of integration and differentiation. 



Lemma 4.2. If for some Eq > and A: G {1, . . . , Af — 1} we have that 
(4.5) \Skiy,z)\<eo, ^y^zeKo, 



then 



\Sk{yr,yr)\<C[j-^j (£o+-B)|lnr|, n = 2or4, 



2fl2iVi 



(4.6) is,(y,,y,)\<c(-^y {so+E)r-\ n = 3, 

^2 o2Ni 

\Sk{yr,yr)\<c(^-^^y {eo+E)r^-", n>5, 

where yr — Pk+i — Ti^iPk+i)- f is small, v{Pk+i) is the exterior unit normal vector 
to dDk at Pk+i and the positive constant C depends on rQ, L, A, B and n. 

Proof. Let the dimension n = 3. We fix z G Kq first and consider v(jj) — Sk{y, z). 
By Lemma [4.11 v solves the equation (—A + q^^^)v = Q \n K f] Wk- Moreover, by 
pTT4|) . we have 

(4.7) |f(y)|<C£; dist(y,Efc+i)-5, y £ K r\Wk. 
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Then, by Proposition 13.91 with 7 = — i, we have, for < 7' < 2ri, 

(4.8) l^'=(2/-'^)l^^(i^j {eo + E)r-i. 
Next, we consider 

(4.9) v{z) = Sk{yr,z), z€KnWk, 

which solves the equation (—A + q^^^)v = in X n Wk, and, by p. 141) . satisfies 

(4.10) \v{z)\<CE {rdist{z,Y.k+i)rK z e K nWk- 

By Proposition [HiH again, we then obtain estimate (|4.6p for n = 3. 

The proof for other dimensions follows from the same proof with a few modifica- 
tions. For 71 = 2,4, a modified version of Proposition 13. 9[ as stated in Remark l3.10i 
needs to be applied. For n > 5, one can apply Proposition 13.91 with 7 = 2—^. D 

Proof of Theorem\2l\ Let 



and 

From the Alessandrini identity (see for instance. Chapter 5 of [9]) 
(4.11) 

{q^^^^q^^^){x)G^{x,y)G2{x,z)Ax = {{Ki-k2)G^{:y),G2{-,z)), yy,zeKo 



and Proposition 13.31 we find that 

(4.12) \Sk-iiy,z)\<C{s + Sk-i). 

Let Pk £ Sfc and yr = Zr ^ Pk — ri'iPk), where v{Pk) is the exterior unit normal 
vector to dDk-i and r is small. We write 

(4.13) Sk-i{yr,yr) = h+h 
with 

(4.14) h= f (g(i)-g(2))(a.)Gi(x,2/,)G2(x,2/.)dx 

JBpo(Pfc)nDfc 

and 

(4.15) h^ f {q^^^-q^^^){x)Gi{x,yr)G2ix,yr)dx, 

JUk-i\iBpg{Pk)nDk) 

where po = ^. 

For 7T, = 3, by Proposition 13.41 we have 

(4.16) Ihl <CE. 
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We estimate /i as follows: 



Gi{x,yr)G2ix,yr)dx 



Spo(Pfc)nDfe 



-^1 (1) (2) I 

>\rk -% \ 



Sp(,(Pfc)nDfc 



r{x,yr)r(x,yr)dx 



(4.17) 



{Gi{x,yr) -T{x,yr))T{x,yr)dx 



Spo(Pfc )n-Dfe 



Bpg{Pk)nDi, 



(G2(x, yr) - T{x, yr))T{x, yr)dx 



(Gi(x,j/,.) -T{x,yr)){G2ix,yr) -T{x,yr))dx 



By Propositions 13.41 and the fact that 



{Gi {x,yr)-r(x,yr))r{x, yr)dx 



Bp^iPk)nDk 

1 
< - 

- 2 



Spo(Pfc)nDfc 



1, 



2\G,ix,yr)-T{x,yr)\' + -|r(x,2/,)|^ ] dx, z = 1,2 



and 



(Giix^yr) ~V{x,yr)){G2{x,yr) -V{x,yr))dx 



Bpg(Pk)nDk 



< 



Bpg(Pk)nDk 



{\Gi{x,yr) -~r{x,yr)\'^ + \G2{x,yr) -r{x,yr)\'^)dx, 



we obtain that 



|/il>l'zi''-gflf ^/ \rix,yr)\'dx~c 



Using the explicit form of r(x, y), we find that 

(4.18) >C|g[')-gf |r-i-C£;. 

Now, by Lemma l¥^ and (|4.12l) . we have 



\Sk-iiyr,yr)\ < C 



£ + Sk-i 



t;p' 



Hence, using (|il^ . (im|) and (|i?T5)) . we have 

£ + 4-1 \ 



(£ + 4-1+^)?^"^ 



\'i'-<i'\r-^<G[ E 



e + 6k-i+E 
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(e + 4-1 + E)r- 



so that 

(4.19) 

Noting that 



imphes 



,i"-a<c(.+..-.+^)((,f±^ 



Tr^ V 7 . Vre (0,2.0 
In (6^ 



""'' > - ^ , , Vre(0,2ri) 



2„2Ni 



</3 



r - 12rihi3' 



we get 

(4.20) |gi'^-gf |<C(£ + 4-i+S) 



e + 4-i +-B 



By taking r 



In 



\e+5t,_i+E) 



-1/4 



and noting that 



/32"i(12riln3)-V2 



(-^-') 

for sonic constant C, we obtain that 

(4.21) |g«-gf|<C(£ + 4-i+i?) 

We let 



< Cr, Vr > 



hi 



£ + 4- 



e + 4-1 + E 



uj{t) = 



lni|-4, G<i<e-3, 
^ 3-i, i>e-3. 

Noting that the function 1 1-> tw„(l/t) is increasing, we have 

£ + 4-1+^ / e + 4-1 



£ + 4-1 V^ + 4-i + ^ 



>c^(l), 



hence 



4-1 <e + 4-i <{Lo{l))-\e + 5k-i+E)u 



g + 4-i 

£ + 4-1 + E 



which with (|4.2ip gives that 



(4.22) 



4 < C(£ + 4-i+£^)w 



£ + 4-1 
£ + 4-1 + E 



The above choice of r is possible only if 

£ + 4-1 



In 



£ + 4-1 + E 
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-1/4 



< 2ri. 



However, if 



In 



e + 4-1 



£ + Sk- 



E 



-1/4 



>2ri, 



that is, 



the fact that 



^ + '^fe-i > e-(2'-i)- 



£ + Sk- 



sup t 

rG(0,2ri) 
t6(e-(2-l)-^l) 



E 



,32"i(12'-lln3)-V^ 



Inih 



is finite shows that (|4.2ip stiU holds true, then (|4.22p follows. 
We iterate (|4.22p . starting from Sq — 0, and find 



(4.23) 



Sk+£<{C + i'''f{E + s) ^u (^) , 



where ujk is the composition of w fc times with itself. We recall that E = 5m , whence, 

e 



(4.24) 
so that 
(4.25) 



E + £<{C + 2,^'''fHE + e) 



COM 



e + E 



- ^^/((C + 3V4)-M) ' 



which completes the proof for dimension n = 3. 

The proof for n = 2 and n = 4 follows from a careful inspection and adaptation 
of the above proof for n = 3. By Proposition 13.51 and 13.41 and the explicit form of 





r{x,y) = 
T{x,y) = 


2tt{x2 - 2/2) 

|a;-y|2 

1 


n = 2 
n = 4 


we obtain that 


4:Tt'^\x — J/ p ' 



(4.26) ki'^-gf |<C(e + 4-i+i?) 



e + 5, 



k-l 



r,.P' 



In) 



Then, by taking r = 



e + 6k-i+E^ 
— rj — g and adapting the function uj{t) according to 



£ + <5fc 



Qit) 



Intl 



< t < e-^, 



t > e- 



we end up with 

(4.27) 



^^ l-^M((g + 2)-*0 ^^^ 



•^J 



((C + 2)-*0 
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which completes the proof for n — 2 and n — A. 

Let us sketch the required modifications of the proof for higher dimensional 
cases(n > 5) below. First, one can use the same decomposition of the singular func- 
tion as in (I4.13[) . (j4.14|) and (j4.15|) . and by Proposition I3.3[ the same upper bound 
estimate of I2 as in (|4.16p is obtained. Then, because the order of the singularity of 
r(-, y) increases as the dimension n increases, G{-, yr) — r(-, j/^) may not be uniformly 
bounded in L'^{Bpg{Pk) fl Dk) with respect to r. A feasible modification here is to 
compare the orders of singularity of G(-, y^) — r(-, y^) and r(-, y^)- More precisely, we 
can estimate 



iBp„{Pk)nDk 
using Holder inequality, as 



{Gi{x,yr) -T{x,yr))T{x,yr)dx 



, 1 = 1,2, 



/. 



{Gi{x, yr) - r(x, 2/r))r(x, yr)dx 



BpgiPk)nDk 

< I \{Gi{x,yr)-T{x,yr))T{x,yr)\dx 

JBpg{Pk)nDt 



Substituting the above inequality into (|4.17|) and noting the positiveness of r(-, j/^), 
we obtain the estimate of the lower bound of /i as 



(4.28) 



\Ii\>\qi'^-q^'^\{m;yr)\\UB„,(P.)nD.) 

- \\G2i-,yr) - r(-, yr)||i2(t^^)||r(-, yr)\\L^{Bp^iPk)nD,) 

~ l|Gi(-,yr) -r(-,y^)||L2(o,) \\G2{-,yr) ~T{-,yr)\\mno) - C). 



By the explicit form of T{x,y) and Proposition 13. 3[ especially (13. 7p . we observe that 
Gi[-,yr) — r(-, j/r) has the lower order of singularity than r(-, yr). Hence, by Young's 
inequality as in the previous proof for n = 3, we conclude that 



|/i|>C|g«-gf|r4-' 



GE 



for r small. Then, following the same argument with the same value of r and noting 
that 



/3=^i(12riln3)-V^ 



holds true for any n > 5, we obtain that 



< Cr 



n—4 



Vr>0 



(4.29) 



ki'^-gf|<C(e + 4-i+i?) 



In 



£ + Sk-i 



^e + Sk-i+E^ 
The last step of the modifications is to adapt the function uj{t) according to 



^it) = 



llntl 



4-T] 

n i 



<t<e- 

t > e-". 
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Then we end up with 

(4.30) E < 



l-^^/((C + n^)-^0 



which completes the proof for n > 5. D 

5. Exponential behavior of the Lipschitz stability constant. In this sec- 
tion, we give a model example to show that the Lipschitz stability constant C = 
C{n,ro, L, A, N) in Theorem 12.71 behaves exponentially with respect to the number 
N of the subdomains. The construction is an analogue of the construction in [11], 
pertaining to the inverse conductivity problem. 

Let n be the unit ball Bi(0) C M" and D = [-1/2, 1/2]" be the cube of side 1 
centered at the origin. We define the class of admissible potentials by 

(5.1) A^ {qe L°°{n) \l/2<q< 3/2 in O and q ^ 1 in n\D} 

and denote the operator from potential q to Ag by F, which maps A into 
C{H^^^{dn), H-'^/'^idQ)). We fix a positive integer N and let iVi be the smallest 
integer such that N < TV". We divide each side of the cube D into A^i equal parts of 
length h = 1/Ni and let Sni be the set of all open cubes of the type 

D' - (-1/2 + {j[ - l)h, -1/2 + j[h) X • • . X (-1/2 + (j-; - l)h, -1/2 + j;», 

where j'l, ■ ■ ■ , j^ are integers belonging to {1, . . . , A^i}. We order such cubes as follows. 
For any two different cubes D' and D" belonging to Sni , we say that D' ~< D" if and 
only if there exists an lo € {1, . . . , n} such that jl = j'- for any i < io and j^'^ < j'/^. 
We define 

N 

AN^{qe L°-{n) I g(x) ^J2'ii^D,{x)+XDo{x), q, £ [1/2,3/2]}. 

Our aim is to estimate from below the Lipschitz constant C{N) in terms of N. A 
simple computation shows polynomial behavior of the lower bound estimate of C{N). 
To obtain the exponential estimate, we then need to employ a topological argument. 
Consider a subset An C An defined by 

^ fl 3 1 

AN = {q£L°°{n) \qix)^Y.'ii^D,ix)+XD„{x), q,£ 3' ^' 2M ' 

It is easy to check that An is a 1/2-net of An with 3^ elements and, for any two 
different qi,q2 £ An^ we have \\qi — q2\\L°°(n) = 1/2- Based on Mandache's result 
[101 Lemma 3], there exist a constant K^ which only depends on dimension n, such 
that for every e G (0,e~^), there is an £-net Y for F{A) with at most e^(~'"^^ 
elements. For e £ (0, e~^) and A^ G N let 

Q(e, AT) = e^^'(- '"")'""'. 
Note that 

3JV>gK(-lne)-- 
20 



if 

where Ki = (K^^ Ina)^/^^""^). There exists Nq such that for iV > A^o we have that 
e < e^^. Thus, for N > No, if we take e = eo we have 3^ > Q{e,N). Then, there 
exist two different qi,q2 £ An such that \\qi — q2\\L=°{n) = 1/2 with their images 
under F in the same bah of radius e centered at a point of Y , that is, 

from which we get 
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